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1. Introduction 



The 0{a?-) initial state QED corrections to e+e annihilation into a virtual photon or Z-boson 
were only calculated once in Ref. [1] in a standard Feynman diagram calculation. Due to the 
smallness of the ratios E, = rn^/s', with m the electron mass and s' the cms momentum squared of 
the virtual gauge boson, power corrections in t, are negUgible and only the logarithmic and constant 
terms in this parameter have to be maintained. The renormalization group equation, applied to this 
process, allows to derive all contributing terms in the on-mass-shell scheme. Our goal is to calculate 
all these terms and to check the result derived in [1] previously. 

In this Umit the differential scattering cross section is given in Mellin space by 
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Here the F's denote the massive operator matrix elements, the a's are massless Wilson coefficients, 
and the P's are splitting functions. a^")(5') is the Born level scattering cross section in terms of 
the invariant mass s of the initial e^e^ pair and L = In [s' jm^). It is convenient to represent the 
differential scattering cross section in terms of three contributions, the flavor non-singlet terms with 
a single fermion Une (I), those with a closed fermion Une (11), and the pure-singlet terms (III). ^ 
These contributions are : 
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The massless Wilson coefficients are known in the literature in the MS scheme for the Drell- 
Yan process [2, 3]. The leading and next-to-leading order splitting functions /^^'^ and p/y ^ are 
well-known from QCD. The only missing piece we need to calculate are the massive fermionic 
operator matrix elements r|*^' to 0(a-^). In the following sections, we describe this calculation. 

^In Ref. [1] four contributions were considered dividing those to I into two pieces according to the genuine 2 — » 3 
particle scattering cross sections. 
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2. Calculation techniques 

To obtain T^})'^ we calculate the diagrams in figure |l] (and the corresponding self-energy 
diagrams) using the rules shown in figure ^ for an external massive fermion (anti-fermion) ^ for 
general integer values of the Mellin-moment n in the form 



(2.1) 



In this way we can separate the splitting function Pa {x) [4] from the remainder non-singlet splitting 
function in the single pole term \/e and in the constant part of the operator matrix element. 



Since we consider the radiative corrections to the time-like process 
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Zo up to two- 



loop order for both initial states, the particle nature has to be conserved, i.e. only e 
e+ — > e+ transitions contribute, which requires a corresponding projection in Y^ee{x)- 



and 



All of the diagrams can be written as a linear combination of integrals of the following type 
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and a couple more of integrals associated with the diagrams that have two photons coming out 
of the vertex. Here D\ = k\ — m^, D2 = k\ — rrp-, D3 = (^1 — pY, D4 = {k\ — ^2)^ and D5 = 
(^2 — ^1 + pY ~ with p the momentum of the electron. For example, the cross-box diagram 
appearing in figure [l] can be written in terms of these integrals as follows 
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These integrals can be calculated in various ways. In case of external gluon lines the 0{aj) 
the integrals appearing in the massive operator matrix elements [6] can be calculated using Mellin- 
Barnes techniques and representations in terms of generalized hypergeometiic functions without 
using the integration-by-parts technique. Working in Mellin space one may solve difference equa- 
tions for harmonic sums [7], cf. [8], which leads to ab initio compact expressions. In the present 



We recalculated the LO OME's and corrected typographical errors in [1]. 

^In space-like processes as deeply inelastic scattering e+ transitions do even occur in 0{a^L?), cf. [5]. 
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calculation we follow the more traditional way of calculation, which was also used in [9]. The mas- 
sive 5-propagator integrals are reduced to 4-propagator integrals using the integration-by-parts 
technique [10]. For example, one of the most complicated integrals coming from the cross-box 



diagram, namely ^^2\ 1 1 > can be decomposed as follows 
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The 4-propagator integrals obey representations in terms of up to three Feynman parameter 
integrals over the unit cube. In course of the calculation we arrange the structure of the Feynman 
parameter integrals such that they have the form 
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(2.8) 



The arts of integration consist now in performing the numerous integrals ( ^ ) in Feynman pa- 
rameter space. This partly requires Mellin-Barnes representations, repeated conformal mapping 
of integration variables, the use of parameter mappings for hypergeometric integrals, and extended 
regularization. Numerous integrals of (Nielsen) poly logarithms [11] normally of complicated argu- 
ments have to be performed. It turns out that the final results can be expressed in terms of Nielsen 
integrals, partly weighted by denominators l/x,l/{l —x^, 1/(1 k <3,l <2. This structure 
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is expected similar to that occurring in case of a wide variety of other 2-loop quantities [12]. For 
example, the integral mentioned above gives 



Detailed checks of the analytic calculations were performed using numerical integrations at high 
precision with MAPLE for fixed moments and individual integrals whenever possible. Mellin- 
Bames integrals were performed for comparison in a series of fixed moments using the code MB 
[13]. We also compared to a few low (0th-2nd) moments for some of the scalar integrals to the 
results obtained using tarcer [14] based on the algorithm [15]. 

We compared the results of our calculation to that in [1] in the logarithmic orders and agree 
after correcting a series of typographical errors in some functions there, decomposing the scattering 
cross section according to the renormalization group. The final result including the constant term 
will be pubhshed soon. 

3. Conclusion 

Since the scale np- /s' <C 1 the QED initial state radiative corrections to e+e^ Z* /y* to 0{a^) 
can be calculated neglecting the power corrections. Alternatively to a direct calculation one may 
use the renormalization group to arrange the calculation in terms of convolutions of splitting func- 
tions, massless Wilson coefficients and process-independent massive operator matrix elements. 
The latter quantities with outer massive fermion lines form the essential part of the present calcula- 
tion. We used the integration-by-parts technique to reduce the integrals to 4-propagator integrals, 
which lead to Mellin moments in terms of triple Feynman parameter integrals. The differential scat- 
tering cross section can be obtained as assembly of the massive operator matrix elements, massless 
Wilson coefficients for the Drell-Yan process and the spUtting functions to two-loop orders. 
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